On the stability of the critical state with inhomogeneous 
temperature in composite superconductors 

N.A. Taylanov 

Theoretical Physics Department, 
Scientific Research Institute for Applied Physics, 
National University of Uzbekistan, 
Tashkent, 700174, Uzbekistan 

E-mail: taylanov@iaph.tkt.uz 
Abstract 

The problem of the thermal and magnetic destruction of the critical state in composite superconduc- 
tors is investigated. The initial distributions of temperature and electromagnetic field are assumed 
to be essentially inhomogeneous. The limit of the thermomagnetic instability in quasi-stationary 
approximation is determined. The obtained integral criterion, unlike the analogous criterion for a 
homogeneous temperature profile, is shown to take into account the influence of any part of the 
superconductor on the threshold for critical-state instability. 
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While dealing with instabilities of the critical state in hard superconductors, the charac- 
ter of the temperature distribution T(x,t) and that of the electromagnetic field E(x,t) are 
of substantial practical interest [1]. This derives from the fact that thermal and magnetic 
distractions of the critical state caused by Joule self-heating are defined by the initial tem- 
perature and electromagnetic field distributions. Hence, the form of the temperature profile 
may noticeably influence the criteria of critical-state stability with respect to jumps in the 
magnetic flux in a superconductor. Earlier (cf., e.g., [2]), in dealing with this problem, it was 
usually assumed that the spatial distribution of temperature and field were either homoge- 
neous or slightly inhomogeneous. However, in reality, physical parameters of superconductors 
may be inhomogeneous along the sample as well as in its cross-sectional plane. Such inho- 
mogeneities can appear due to different physical reasons. First, the vortex structure pinning 
can be inhomogeneous due to the existence of weak bonds in the superconductor. Second, 
inhomogenety of the properties may be caused by their dependence on the magnetic field H . 
Indeed, the field H influences many physical quantities, such as the critical current density 
j c , the differential conductivity ad, and the heat conductivity k. 

In the present paper, the temperature distribution in the critical state of composite su- 
perconductors is investigated in the quasi-stationary approximation. It is shown that the 
temperature profile can be essentially inhomogeneous, which affects the conditions of initi- 
ation of magnetic flux jumps. 

The evolution of thermal (T) and electromagnetic (E, H) perturbations in superconductors 
is described by a nonlinear heat conduction equation [3], 
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v^ = V[KVT]+jE, (1) 

a system of Maxwell's equations, 

-* ldH 

r0tE =-clf (2) 
rotH = —j (3) 

and a critical-state equation 

j=j c (T : H)+j r (E). (4) 
Here v = v(T) is the specific heat, k = k(T) is the thermal conductivity respectively; j c is 

— * 

the critical current density and j r is the active current density. 

We use the Bean-London critical state model to describe the j c (T, H) dependence, ac- 
cording to which j c (T) = jo — a(T — T ) [4], where the parameter a characterizes thermally 
activated weakening of Abrikosov vortex pinning on crystal lattice defects, jo is the equilib- 
rium current density, and To is the temperature of the superconductor. 

The dependence j r {E) in the region of sufficiently strong electric fields (E > Ef \ where Ej 
is the limit of the linear region of the current-voltage characteristic of the sample [2]) can 

TIC? <7 n H c 

be approximated by a piecewise-linear function j r m (JfE, where cry = ~ "^ C2 is the 

ivhc 

effective conductivity in the flux flow regime and n is the viscous coefficient, $o — is the 
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magnetic flux quantum, a n is the conductivity in the normal state, H C2 is the upper critical 
magnetic field. In the region of the weak fields (E < Ef), the function j r (E) is nonlinear. 
This nonlinearity is associated with thermally activated creep of the magnetic flux [5]. 

— * 

Let us consider a superconducting sample placed into an external magnetic field H = 

dH 

(0, 0, H e ) increasing at a constant rate = i/=const. According to the Maxwell equation 

at 

(2), a vortex electric field E = (0, E e , 0) is present. Here H e is the magnitude of the external 
magnetic field and E e is the magnitude of the back-ground electric field. In accordance with 
the concept of the critical state, the current density and the electric field must be parallel: 
E || j; where H e is the amplitude of the external magnetic field and E e is the amplitude of 
the external electric field. 
The thermal and electromagnetic boundary conditions for the Eqs. (l)-(4) have the form 

dT 



dx 



dE 

dx 



+ w [T(0)-T } = 0, T(L) = T , 



=o 



(5) 



0, E(L) = 0, 



x=0 



For the plane geometry (Fig.l.) and for the boundary conditions H(0) = H e , H(L) = 0, 

cH e 

the magnetic field distribution is H(x) = H e (L—x), where L = - — - is the depth of magnetic 
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flux penetration into the sample and wq is the coefficient of heat transfer to the cooler at 
the equilibrium temperature T . 

The condition of applicability of Eqs. (l)-(4) to the description of the dynamics of evolution 
of thermomagnetic perturbations are discussed at length in [1]. 

In the quasi-stationary approximation, terms with time derivatives can be neglected in 
Eqs. (l)-(4). This means that the heat transfer from the sample surface compensates the 
energy dissipation arising in the viscous flow of magnetic flux in the medium with an effective 
conductivity 07. In this approximation, the solution to Eq. (2) has the form 



E — —(L — x). 

c 



(6) 



Upon substituting this expression into Eq. (1) we get an inhomogeneous equation for the 
temperature distribution T(x,t), 

d 2 Q 



dp 2 



P© = f(p). 



(7) 



Here we introduced the following dimensionless variables 



f{p) = -[l + rup] 



Jo 

aT ' 



= 



T-T 
To 



P = 



L — x 



and the dimensionless parameters uo = , and, r = ( — : \ , where r characterizes 

cjo \aHL 2 / 
the spatial scale of the temperature profile inhomogeneity in the sample. Solutions to Eq. 

(7) are Airy functions, which can be expressed through Bessel functions of the order 1/3 [6] 
0(p) = C lP ^K 1/s (^p 3 / 2 ) + C 2 p^I 1/3 (Jp 3 / 2 ) + o (p), (8) 

0o(p) = p^K 1/3 (Jp 3 / 2 ) j\l + r^pTh/s (Pr) d Pl - 



P 1/2 h/s (Ip 3/2 ) + (jpf ) d Pl , 

where C\ and C2 are integration constants, which are determined by the boundary conditions 
to be 



-w L&(0) + ^ 
dp 



Ci = 0, c 2 



p= — 

r 



L\ 1 ^ 2 /2 

W \ - h/3 ^P~ 3/2 



From the Maxwell equation (2), the temperature inhomogeneity parameter can be expressed 
in the form 
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a 



aH 2 Ht K 



1/3 



(9) 



It is evident that a ~ 1 near the threshold for a flux jump, when 



aH]_ 



1, even under the 



quasi-stationary heating condition -— - << 1; where t K = is the characteristic time of 



K 



the heat conduction problem. 
Let us estimate the maximum heating temperature O m in the isothermal case w — — > 1. 

Lj 

The solution to Eq. (7) can be represented in the form 

{x %rn) 



Q(x) = <S>m - PO- 



near the point at which the temperature is a maximum, x = x m (Fig.l). 
With solution (10) being approximated near the point x m = — with the help of the 



thermal boundary conditions, the coefficient po can be easily determined to be 
and the temperature can be written as 



LP- 



e(x) = e r 



4 / L 



(ii) 



Substituting this solution into Eq.(7), the superconductor maximum heating temperature 
due to magnetic flux jumps can be estimated as 



Jo 



a f H 



(L - x m ) 



H 
1 

ckT 



[L - x, 



L 2 CK 



(12) 



For a typical situation, when — << (L — x m ) the estimation for m is 



jo + ^-{L - x m ) 



HL 2 

ckT 



(L — x m ). 



(13) 



Here, the parameter 7 ~ 1 (for a parabolic temperature profile 7 ~ 8). It is easy to verify 
that for typical values of jo = 10 6 A/ cm 2 , H = 10 4 G/s, and L = 0,01 cm the heating is 
sufficiently low: O m << 1. In the case of poor sample cooling, w — 1 -j- 10erg/(cm 2 sK), the 



9 m is 



Hj L 2 



0,5; 



CWqTq 

i.e., the heating temperature can be as high as ST m = ToO m ~ IK. One can see that 
in the case of poor sample cooling, the heating can be rather noticeable and influences the 
conditions of the thermomagnetic instability of the critical state in the superconductor. 
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Let us investigate the stability of the critical state with respect to small thermal 5T 
and electromagnetic 5E fluctuations in the quasi-stationary approximation. We represent 
solutions to Eqs. (l)-(4) in the form 



T(x,t) =T(x) + exp\ — 



At 

T K 



5TlZ 



(14) 



E(x,t) = E(x) + exp ^ — 



At 

T K 



where T(x) and E(x) are solutions to the unperturbed equations obtained in the quasi- 
stationary approximation describing the background distributions of temperature and elec- 
tric field in the sample and A is a parameter to be determined. The instability region 
is determined by the condition that ReX > 0. From solution (14), one can see that the 

characteristic time of thermal and electromagnetic perturbations tj is of the order of t K /X. 

5T 5E A , . , c „ . 

we obtain the following 



Linearizing Eqs. (l)-(4) for small perturbations y 
equations in the quasi-stationary approximation: 



T(xY E(x) 



« 1 



u = T?1^ + [j{x) + afE{x)] 5E ~ aE{x)6T ' 



^^ Wf SE-aS T] 



(15) 



L 2 dx 2 c 2 t 



Eliminating the variable 5T between Eqs. (15), we obtain a fourth- order differential 
equation with variable coefficients for the electromagnetic field 5E: 



d 4 5E 
dz 4 



A(l + r) + 



E(z) 



d 2 5E 
dz 2 



+ A[Ar - B(z)\5E = 0. 



(16) 



Here, we introduced the following dimensionless variables: 



x 



z = 



AiraL 2 
c 2 u 

K 



v B{z) - 



aL 2 ' 



j(z), j(z) = a f E(z)-[ Jo -a(T(z)-T )], E(z) 



HL 



1-z), 



crz/ 



E K 



ToJ 



One should keep in mind that the variable T(z) and E(z) are given by Eq. (8), in which 

p = — (1 — z). Using the relation between 5T and SE given by Eqs. (15), we write the 
r 

boundary conditions to Eq. (16) in the form: 



d 2 5E 
~dz l2 ~ 



0, 



6 El 



z=l 



z=l 



o, 



d 3 5E 
~dz^ 



= -W 



d5E 

dz 



z=0 



0. 



'd 2 5E 
~dz" 2 ~ 



\t5E 



2 = 



(17) 



z=0 



■5 



where W = is the dimensionless thermal impedance. 

K 

The condition for the existence of a nontrivial solution to Eq. (16) subject to boundary 
conditions (17) allows one to determine the boundary of the critical-state thermomagnetic 
instability in a superconducting sample. This problem is complicated, and its analytical solu- 
tion cannot be found in a closed form. We will consider the development of thermomagnetic 
instability in the dynamical approximation, which is valid for composite superconductors 
with high value of <JfE. 

The dynamical character of the instability development leads to the predominance of heat 

D K 
diffusion over magnetic flux diffusion in the sample: r = — — >> 1 [1], where D t = — and 

C 2 

D m = are the coefficients of the thermal and magnetic diffusion, respectively. 

In this case, as seen from Eq. (14), the characteristic times tj of temperature and electro- 
magnetic field perturbations have to satisfy the inequalities tj » t K (A << 1) and tj «t m 

L 2 L 2 

(Ar >> 1), where t K = — and t m = — — are the characteristic times of the thermal and 

E>t Dm 
magnetic diffusion, respectively. 

As well known that [1], the effective magnitude of conductivity 07 is greater in composite 
superconductors than the one in hard superconductors. We can assume that induced normal 
current (jfE compensates the decreasing of critical current j c (T), caused by increasing of 
temperature and obviously prevents magnetic flux penetration into the sample. In this case 
we can neglect the moving of magnetic flux. In the other word, thermomgnetic instability 
developes slowly as compared with the heat diffusion with characteristic time of increasing 
tj ~ % » t K (A « 1) or t m » ^ = % = tj (Ar » 1). 

A AT A 

In the approximation (r >> 1, Ar >> 1, A << 1) Eq. (16) is reduced to a lower order 
differential equation 

d 2 5E 



dz 2 



+ 



5E = 0. (18) 



In the case of r >> 1, the instability threshold depends on the electrodynamic boundary 
conditions at the surface of the sample only slightly. Therefore, the electrodynamic boundary 
conditions at the boundaries of the current-carrying layer (z — 0, z — 1) can be neglected 
and one can keep only the thermal boundary conditions to Eq. (18). 

Multiplying Eq. (18) by 5E and integrating the result with respect to z over the interval 
< z < 1, we obtain 

1 1 
B(z)5E 2 dz - 




1 

A = ^ j u — : (19) 



5E 2 dz 



where we use the equality 
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and the boundary conditions. The right-hand side of Eq. (19) has a minimum at A = A c : 

1 } 

- / B(z)5E 2 dz 

Ac = ^4 , (20) 

5E 2 dz 



o 



then 

} (d5E\ 2 



1 ? j \ dz ) 
- / fl(*)n|d* = ^ . (21) 



Here we introduced the following unit vector 



SE 2 
5E 2 dz 



n \ = ^ 1 , Jn 2 E dz=l. 



Since we do not know the function SE(z), we try, following [7], to obtain an integral 

estimation of the instability growth increment and the low boundary of its occurrence. The 

H L 

behavior of the integrand in Eq. (21) is basically determined by the factor E = (1 — z), 

which is equal to zero at z — 1 (the other factors change more smoothly). Hence, the 
integrand reaches its maximum at z = and the upper estimate for A c is 

A C <^. (22) 

It is evident that A c << 1 and A c r » 1 at r » 1. Numerical evaluation gives A c ~ 
0.1 -=- 10~ 2 at t = 10 3 . 

Equations (21) and (22) enable one to write the instability occurrence criterion in the form 

1 ) J [dz 

- J B{z)n 2 E dz > ^ . (23) 



SE 2 dz 

This criterion essentially depends on the boundary conditions and the functions j(z), E(z), 
T{z). Figure 2 presents graph of the function T(z). Inequality (23) can be strengthened by 
means of an evaluation, 
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v 7 

which can be easily verified by expanding the function SE(z) in a Fourier series: 

, n/ > , 7vz(2m + 1) 
5£(z) = A m cos — y — '- . 

Let us now try to strengthen inequality (23) further. For this purpose, we consider the 
integral 

1 21 

J 9{z){n 2 E - l)dz = J g{z){n 2 E - l)dz + J g{z){n 2 E - l)dz. 

zi 

The last term can be represented in the form 

i i 
/ g(z)(n 2 E - l)dz = (g + - g_) j g(z)(n 2 E - l)dz, 



taking intermediate values of the g- in the range z < z\ and g > g- in the range z\ < z < 1 
outside the integral. It is evident (Fig. 3) that 

i 

g(z)(n 2 E -l)dz<0 



o 

or 



J g{z)n 2 E dz < j g(z)dz. (25) 
o o 

With inequality (25), the instability occurrence criterion can be represented in the form 



/ 



7T 2 



B{z)dz>—T. (26) 
o 

Inequality (26), unlike the analogous criterion for a homogeneous temperature profile, has 
an integral character and takes into account the influence of each part of the superconductor 
on the threshold for the superconducting-state instability. If condition (26) is satisfied, 
then small fluctuations of temperature ST and electric field 5E in the superconductor will 
exponentially increase with time. The most probable result of the development of such an 
instability would be a transition from a critical state to a resistive one. 

We should emphasis that this result was obtained for an arbitrary temperature dependence 
of thermophysical parameters v and k of superconducting material and for an arbitrary 
function j(H). Moreover, since the system of equations (l)-(4) is invariant with respect 
to an arbitrary translation, the wave propagation condition can be found for an arbitrary 
critical current density dependence on T and H. 
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FIGURE 

Fig.l. The distrubation of the temperature profile Q(x). 
Fig.2. Plots of the function T(z). 
Fig. 3. Plots of the function g(z). 
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